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We compute the three loop (3 function of the Wess-Zumino model to motivate implicit regularization 
(IR) as a consistent and practical momentum-space framework to study supersymmetric quantum 
field theories. In this framework which works essentially in the physical dimension of the theory 
we show that ultraviolet are clearly disantangled from infrared divergences. We obtain consistent 
results which motivate the method as a good choice to study supersymmetry anomalies in quantum 
field theories. 
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I. INTRODUCTION 

Dimensional regularization (DR) is a remarkable 
framework which besides preserving gauge invariance is 
relatively simple from the computational standpoint. For 
this reason it has become the stantard method in pertur- 
bative calculations in quantum field theory. The regu- 
larization and renormalization of supersymmetric gauge 
theories is, however, a more involved problem. Mo- 
mentum space subtraction schemes such as BPHZ 
although supersymmetric are not gauge invariant 
DR does not mantain the balance between bosonic and 
fermionic degrees of freedom due to the analytic contin- 
uation on the space-time dimension and therefore breaks 
supersymmetry. Such breaking demands the calculation 
of compensating supersymmetry restoring counterterms. 
A practical useful modification of DR is dimensional re- 
duction (DRed) 3] which is however mathematically in- 
consistent and cannot work at all orders. As long as 
such inconsistencies can be tamed and symmetry restor- 
ing counterterms can be unambiguously generated by im- 
posingthe validity of the Slavnov- Taylor and Ward iden- 
tities |j] the use of DR and DRed are obviously justified. 
However it is not always the case that the symmetry of 
the Lagrangian is still a symmetry of the full quantum ef- 
fective action. Supersymmetry anomalies can in principle 
be generated and some erroneous claims about their ex- 
istence have occurred because it is difficult to distinguish 
between a genuine anomaly and an apparent violation of 
a supersymmetric Ward identity due to the use of an ill- 
defined regularization scheme [5| . In the case of a genuine 
anomaly no symmetry restoring (finite) counterterms can 
be obtained. 

A regularization/renormalization framework that: a) 
does not modify the field theoretical content of the bare 
Lagrangian (and hence does not unnecessarily complicate 
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the Feynman rules) ; b) works in the physical dimension of 
the model; c) preserves gauge invariance and d) is friendly 
from the calculational viewpoint, is therefore desirable to 
tackle the issues that we discussed above. In this sense 
two relatively new frameworks deserve special attention: 
Differential and Implicit regularization/rcnormatization 
(DfR and IR, respectively). 

DfR is a coordinate space method that defines the 
correlation functions without introducing a regulator or 
counterterms j(aj-jl6i]. A well defined prescription that (in 
a minimal sense) extends product of distributions into 
a distribution automatically delivers renormalized finite 
amplitudes. The latter contains an arbitrary mass scale 
which must be introduced by dimensional reasons and 
plays the role of a renormalization group scale. Gauge 
invariance may be sistematicaly implemented in a con- 
strained version of DfR, at least to one loop order. Fi- 
nally contact with momentum space is made by means 
of Fourier transforms. 

On the other hand IR is a momentum space scheme 
where scattering amplitudes with fixed external momenta 
are constructed in the first place. In section [H] we briefly 
outline this method. 

The purpose of this paper is to motivate IR as a sound, 
symmetry preserving and practical framework to study 
dimension specific theories, among which chiral, topo- 
logical and supersymmetric gauge theories are of prime 
interest. We use the massless Wess-Zumino model as a 
testing ground for the consistency of IR in preserving 
supersymmetry. We present the Feynman diagram cal- 
culation of this model in sections HVl IVl and lVIl As a non- 
trivial check we compute in section lVlIl the /3-function to 
three loop order and verify the agreement with other con- 
sistent methods. As a by-product we show how infrared 
and ultraviolet divergences are clearly disantangled (in 
opposite to dimensional methods). This is important in 
the study of certain supersymmetry anomalies. We also 
define what is meant by a minimal subtraction within 
IR and compare with DR and DfR. Finally, in section 



IVIIII we address some theoretical and phenomenological 
problems where IR could be useful. 
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II. IMPLICIT 
REGUL ARIZ ATION /RENORM ALIZ ATION 

The main idea behind IR is to isolate the divergences 
from an amphtude as irreducible loop integrals (ILI) 
which do not depend upon the external momenta by ju- 
diciously using the identity: 



1 



[(fc + fci)2 - m2 



N 



{-iyikf + 2h-ky 



(1) 

In the equation above, ki are the external momenta and 
iV is chosen so that the last term is finite under integra- 
tion over k. For instance take the logarithmically diver- 
gent integral in four dimensions 



1 



r)4 (fc2-m2)[(fc+p)2-m2] 
We use ^ with = 1 to write 



(2) 



df^k 



1 



(27r)4 (fc2 - m2)2 
d^k + 2p-k 

(27r)4 (fc2_TO2)2[(fc+p)2-m2 



(3) 



Note that the second term on the rhs of the equation 
above is finite whereas the first is an ILI, namely. 



1 



(27r)4 (fc2 _ ^2)2 



(4) 



which characterizes the ultraviolet behaviour of the am- 
plitude and need not to be explicitly evaluated. They can 
be fully absorbed in the definition of the renormalization 
constants. In order to define a mass independent scheme 
(or in the case of massless theories) we use the identity 
(|18|l which we introduce in section Hvl Such identity in- 
troduces naturally an arbitrary scale which plays the role 
of renormalization group scale. Infinities of higher order 
are equally displayed as ILI such as Iquadim?), Innirn^) 
etc., for quadratic and linearly divergent integrals, re- 
spectively. Local arbitrary counterterms can show up in 
IR as (finite) differences between divergent integrals of 
the same superficial degree of divergence which can sis- 
tematically be cast into a set of "consistency relations" 
according to the space-time dimension. At the one loop 
level they are related to momentum routing invariance in 
the loop of a Feynman diagram. A constrained version 
of IR corresponds to setting such consistency relations to 
vanish. In this way, abelian and nonabelian gauge invari- 
ance can be shown to automatically implemented. Ulti- 
mately, in a more general context, such local arbitrary 
counterterms parametrized by the consistency relations 
should be left arbitrary till the end of the calculation 
when the symmetry content of the underlying model may 
fix its value. We address the reader to ref. for 
details and applications of IR. 



Overlapping divergences which are the chief complica- 
tion in any renormalization scheme can also be handled in 
a schematic way within IR. This was illustrated to n-loop 
order in the context of <j)'^ theory in 6 dimensions. 



III. PERTURBATION EXPANSION OF THE 
WESS-ZUMINO MODEL 

We closely follow the notation and the conventions of 
ref. [23. The Wess-Zumino (WZ) model superspace ac- 
tion reads 



S 




(5) 



where {x" , 9" , 9") , a = 1,2,3,4, a = +-, a = +,_- are 
coordinates of d = 4, iV = 1 superspace and (p {(p) is a 
chiral (antichiral) superfield. 

The Lagrangian can be written in terms of component 
fields 

cPo - exp [ie^e^ da) X 

X (l/V2{A + iB) + 9"iPa + 1/V2{F - iG)9^) , (6) 

which after eliminating the auxiliary fields F, G, will 
have the form 



C 



1 



1 



.90 



1 



--07 • dip 



gt,{A' + B'Y + ^^{A + i^,B)i;. (7) 



16 



The WZ model involves only one renormalization con- 
stant since only propagator-type diagrams can diverge. 
To see that, one may eraploy the superfield power count- 
ing rules described in [2g to conclude that the three- 
point function is finite. Defining the renormalization con- 
stants as 



3o 



Zg9 



Z 



1/2, 



(8) 



and using that 50^0 finite (and hence ZgZ'^J'^ = 1) 
enables us to write 



'7-3/2 



(9) 



Thus for computing the /3-function it is sufficient to cal- 
culate Zcj, (i.e. the divergent structure of the two-point 
function). We expand into a power series in the cou- 
pling constant. 



(47r) 



Zi 



47r 



47r 



(10) 
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and write the part of the effective action which is hnear 
in (h and (h hke 



_$i~k,e,mk,0,9)A-\k), (11) 



where and Jg stand for J d'^k/(2n)'^ and / (PO^ re- 
spectively. Following we notice that in the renor- 
malized Lagrangian appears only in the kinetic term 
4'o4'o = which amounts to introducing a factor of 

Z'^^ for each propagator in the diagrams. Hence we can 
write 



2 V47r 

6 ,-8 



2 V47r 



(12) 



in which / {A-k)^ F [p^] represents the contribution of 
the one loop diagram (figure ^1, ^g^ / {Att)'^G[p'^) repre- 
sents the two loop contribution (figurelSJ, 5 being a sym- 
metry factor, and g / [Att)^ H ijp') refers to the three loop 
diagrams (figures |21 to EJ). Now we take equation pU|) 
in equation p2|) and reorganize the power expansion to 
write 



= -11 (ty{-p,6,e)(i3{p,9,e){Z^ + 

(13) 




The Feynman rules for the Wess-Zumino model are 
well-known (please see |2^, j23|) so we shall not derive 
them here. The physical interpretation of equation (|13|l is 
straightforward. It defines the renormalization constants 
ZiS at each loop order after the infinities corresponding 
to subgraphs of previous orders are duly subtracted. Di- 
agrammatically it amounts to substituting the divergent 
subdiagrams of a diagram with their finite part order by 
order, as we shall sec in the next sections. We shall per- 
form such substitution after subtracting the divergences 
in a minimal sense within our approach, namely by sub- 
tracting irreducible loop integrals. 




FIG. 1: Fij? 




FIG. 2: G(p2 




FIG. 3: Hi(j? 



IV. ONE LOOP CONTRIBUTION 

The one loop contribution to the propagator correc- 
tion is represented by the diagram depicted in figure ^ 
Application of the Feynman rules to this diagram yields 
25], 26]: 



r^'^ = {-^?9' I d^9^i-p,0)^{p,e) 



1 



From the equation above and H13|l we identify 
F(p2) -i(47r)2 ^ ^ 



(14) 
(15) 



in which we have introduced an infrared cutoff fi. Follow- 
ing our approach, an irreducible loop integral can be sep- 
arated from the amplitude above with the help of equa- 
tion (Q) to give 



F(/) 




(fc2-M2)2[(p+fc)2-^2] 



(16) 



where, henceforth we define 



b = 



(47r)2 



(17) 



and hogi^J''^) is given by equation I^J. Since the limit 
where /i ^ is infrared ill-defined in Iiogdi^), the correct 
ultraviolet behaviour is obtained by exchanging the latter 
against Iiog{X^) provided 7^ by means of the identity 



Ilog{fl^)=Ilog{X^)+bln(^^^ . 



(18) 
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FIG. 4: Hiip'^ 



FIG. 6: He{p'^ 




FIG. 5: Hsip") 



Hence (|18f) splits the ultraviolet and infrared divergences 
and as a byproduct it parametrizes the arbitrariness in 
separating the divergent from the finite part. That is 
because the infrared divergent piece in the right hand side 
of H18I) cancels out against another piece coming from the 
(ultraviolet) finite part of the amplitude which may be 
easily integrated to give 



(19) 



Consequently A is the natural candidate for a renormal- 
ization scale from which we may construct a Callan- 
Symanzik renormalization group equation. A minimal 
subtraction within our approach amounts to defining, 
with the help of l(T^ . 



Zi 



-(47r)2/,„g(A2), AVO, 



(20) 



i.e., we have subtracted only the irreducible loop integral. 
For future reference, we define the finite part of F(j>^) as 



Ffin = -i(47r)261n 



(21) 



where A plays the role of an arbitrary local counterterm. 
Moreover satisfies a Callan-Symanzik renormaliza- 
tion group equation with scale A |23| . 

For a massive theory a minimal mass independent 
scheme is defined in a similar fashion. In |2ll | we have 
compared our minimal subtraction scheme with the MS 
scheme in dimensional renormalization as well as in dif- 
ferential renormalization. The arbitrary scales appearing 
in each framework are related to each other as we shall 
discuss in section IVIll fsee also H^)- 



to (j) (there is no (jjcj) or propagators). The correspond- 
ing amplitude can be constructed by direct application of 
the Feynman rules. After some straightforward algebra 
it can be reduced to 



jp 

^'^^Lk fcV(fc + 9)2(p + g)2 ' 

which enables us to identify 

1 r 1 



(22) 



G{p') = (4^)^ 



(23) 



Let us have a closer look at the O(g^) term in equation 
(P|l viz.. 



(47r) 



4 r 



1 



G(p2) - z,F{p^ 



(24) 



which graphically corresponds to figure |21 Perturbative 
renormalization is inductive in the definition of the coun- 
terterms. Thus in defining Z2 we firstly ought to subtract 
the subdivergence through Zi which was defined in the 
previous order as prescribed in H24|l . We can explicitly 
evaluate l/2G{p^) - ZiF{p^) using ^ and ^ 

to conclude that such operation amounts to substituting 
the subintegration over k in (ESJl with ffi„(p2) divided by 
— i(47r)^ to define: 



G{p'') = biAnY 



1 



q'^{p+q)^ 



In 



-G(p2) - Z.Fip^ 



(25) 



Hereforth the tilde means that we have subtracted all 
the subdivergencies. The graphical interpretation of such 
procedure is clear. 

We proceed to define the renormalization constant Z2. 
Introducing an infrared cutoff in (|25(l enables us to write 



V. TWO LOOP CONTRIBUTION 

The propagator receives just one two-loop contribution 
(figure 121) because the only propagator is the one from (j) 



G(p2) = 6(47r)^ 



(g2 -/i2)[(p + g)2 -^2- 



In 



2^2 



-A^e 



(26) 

The irreducible loop integral can be separated just as we 
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did at the one loop level as follows 
G(p2) = 6(47r)^ ^ ^ 



In 



X In 



2p-q 



9 - Ai' 



(27) 



Notice that the second term on the right hand side is 
finite whilst the first term is an ILI, namely 



1 



(g2 _ ,,2-12 



In 



I'dM)- (28) 



Again, a bona fide renormalization constant should be 
infrared finite. Thus we subtract with A 7^ by 

parametrizing the infrared divergence by means of the 
identity 



(2), 



1 1 2 / M 



- 2 



(29) 



Relation H29|l is the two loop analog of 118(1 . Now taking 
(|29|l and (|28|l into l|27|) permits us to subtract the diver- 
gence by defining the renormalization constant of order 



5" as 



^2 = -^(4.)^/g(A^). 



(30) 



In the next section we shall evaluate the three-loop con- 
tributions to the propagator and define the corresponding 
counterterms. Thus it is convenient to simplify the finite 
part of the two loop contribution as there will be two loop 
subgraphs at three-loop order. We show in appendix A 
that 



Gfi„(p2) = Gfin(p') 

■In 



1. 2/-A^"' 



-In 



-A^e^ 



(31) 



iJa . . . iJg being the result of integrating over the 8 func- 
tions of the superspace coordinates and eliminating the 
covariant derivatives. The factor — i/(47r)^ appears so 
that the iJ^'s here agree with the definition expressed in 
equation H13() . In the latter iJ stands for Hj, + . . . + iJg. 

Take the first three loop diagram (figure . It con- 
tributes to 1(32(1 with 



i(47r)« 



1 



1 



g g2(p + q)2 k-^(k^v + q? 



(33) 



Note that it contains two one- loop subdiagrams. Accord- 
ing to what we discussed earlier we can substitute the 
integral over / with ^fin((7^)/[— i(47r)^] (equation (|^ ) at 
once. As for the integral over k , it can be cast with the 
help of ^ as 




k^{k + q)^{p + k + q)^ 



(34) 



This is useful because doing so we have contributed to 
free the ultraviolet divergent piece from external momen- 
tum dependence. Note that only the first piece of equa- 
tion namely 



i(47r)^ 
8 



1 



1 



1^{p + 9)^ Jk k^ik + q)^ 



bin 



(35) 



is ultraviolet divergent. Thus we may promptly substi- 
tute the integral over k (which represents a one-loop 
subdiagram with external momentum q) in 1(35(1 with 
Ffi„(q2)/[_i(47r)2] to define: 



i(47r) 



q^(jp + qY 



bin 



(36) 



It is worthwhile mentioning that Gfin(p^) is infrared safe 
as it should, since the limit where ^ — )■ is well defined 
through a cancellation of terms in equations l(27() and 



VI. THREE LOOP CONTRIBUTIONS 

The diagrams depicted in figures 3 to 6 represent the 
three loop order contributions to the propagator. Feyn- 
man rules can be directly applied to give [25j : 



.(3) 



(4^) 



^ [n^ip') + H,ip^) + H,ip') + Heip')) ,(32) 



To define the irreducible loop integral which will con- 
tribute to we introduce an infrared cutoff and re- 
move the external momentum dependence in H^{p^) in a 
similar fashion to equations ((26(1 - 1(27(1 . That is to say 



i(47r)^ 



^3(p',m' 
i(47r)2 



(g2 _ ^2)2 



4'i(A*')+•^3(p^M^), 



62 In^ 



.-AVx 
Kq^-fi-^J 



(37) 



in which we defined another logarithmically divergent 
ILI, 



1 



(g2 _ ^2)2 



ln2( 



q - ^J■■ 



r) , (38) 



6 



whereas ^^{p^jfi^) stands for an ultravfolet finite piece. 
This diagram will contribute to with an infrared finite 
term 



.(3) 



i(47r)2 



A 7^ since one can easily verify that 



8 1,3 
- + - In'^ 
3 3 



(39) 



(40) 



In analogy with our calculations at the one and two loop 
orders, the infrared divergent piece in equation above is 
expected to cancel the infrared divergence in TzijP' , 
to render a well defined finite part. This is indeed the 
case as one can prove after some straightforward algebra. 
Therefore we may write 



i?3(/) 



i(47r) 



i[f>,{\^) + g^{p\\^). (41) 



The second three-loop contribution is represented by 
figure 01 from which one gets 



Hi 



i(47r6) 



1 



1 



kkHk + iy 



(42) 



One can see from this graph that it contains clS Oil sub- 
diagram the two loop graph shown in figure |21 which is 
represented by the integrals over I and k in the ampli- 
tude above, q playing the role of external momentum. 
The procedure is identical as we did for H3 so we shall 
only summarize the steps. We replace the integrals over 
I and k by Gfin(Q^)/(47r)^ (equation (|S3) and introduce 
an infrared cutoff, fi. Next we expand the propagator 
which contains the external momentum p in the usual 
fashion so to define an irreducible loop integral which is 
p-independent. We obtain 



H^ip') 



1 



-i(47r)^ 



In 



+ Ti{p^^i'), 



(43) 



being an ultraviolet finite term. Using , , ll^ 
and H4UII allows us to write 



Hi{p^ 



'i(47r) 



■(^/S(A^) 



+ /S(A^) 



(44) 



where Q^lp^^X^), \^ 7^ 0, is now both ultraviolet and 
infrared finite . Finally this diagram contributes to 
with 

Z^ = ^(^I^i^)^I^i^)). (45) 



The diagram displayed in figure [S] is easy to evaluate. 
It also contains two one-loop subdiagrams just as the 
diagram in figure 01 Therefore we expect that it shall 
contribute to Z3 with a term proportional to -^/^^(A^) as 
well. For the sake of completeness we write its contribu- 
tion to (jS21l: 



1 



1 



4 .lqq^{p + qf Jil^{q + l? 
1 



(46) 



Notice that the integrals over / and k in the equation 
above represent the one loop subdiagrams. We proceed 
in a similar fashion as we did for the diagram in figure |3I 
to obtain 

H,{p') = ^^/S(A^) + Q,{P\ X') , (47) 

where the notation is now obvious. Hence, the renormal- 
ization constant Z3 should contain the following contri- 
bution, 



^f-^/S(A^). 



(48) 



in order to cancel the divergence that stemmed from this 
diagram. 

Finally we study the diagram in figure [HI In oppo- 
site to the previous three loop diagrams, it contains no 
subdiagrams at all. Its contribution can be cast as [2fil |: 



i(47r) 



1 



2 Ji.,.kq'np-ini-qni-k)Hq-kr- 

(49) 

An irreducible loop integral can be displayed after ex- 
panding the propagator that contains p in Hq according 
to CP) to obtain 



Ha 



i(47r) 



1 



2 .ii,,,kq'kHm-qni-kriq-kr 

2\ — zjoo 



W) = H^+W) 



(50) 



J-q{p^) being a ultraviolet finite piece as usual. We show 
in appendix B that 

Hq = -il{^Tlfai)Ilog{\^) + e6(/, A^) , (51) 

where Ci^) stands for the Riemann zeta-function. There- 
fore we define the last contribution to Z3, 



Zf = 3z(47r)«C(3)/;o3(A^) 



to finally write 



Z3 - Z3 



(3) 



z. 



(6) 



(52) 



(53) 
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VII. THE /3-FUNCTION 

In pH we compared implicit, dimensional, differen- 
tial and BPHZ renormalization and shown, within these 
intrinsically distinct frameworks, how renormalizations 
schemes and scales are related. It goes as follows. Mini- 
mally subtracting the infinities in dimensional regulariza- 
tion (i.e. removing only the poles) delivers a finite (non- 
counterterm) term which depends upon the renormaliza- 
tion scale fi. It results that such term is the same that 
appears should we employ differential renormalization, 
except for a simple rescaling of its typical arbitrary scale 
(that is to say, a finite counterterm). The latter is the 
scale of a Callan-Symanzik renormalization group equa- 
tion satisfied by the vertex function Q . A minimal renor- 
malization scheme within implicit regularization corre- 
sponds to subtract the irreducible loop integrals in a mass 
independent fashion through relations like IjlSI) , H29I) and 
(|40|l where an arbitrary mass scale A is introduced. It 
turns out that the resulting finite part can be identified 
with the one from differential renormalization (Fourier 
transformed to momentum space) after a simple rescaling 
of A. The difference is that in our framework we do have 
counterterms from which we may calculate the renormal- 
ization group functions. Moreover our renormalizcd am- 
plitude also satisfies a Callan-Symanzik renormalization 
group equation governed by the arbitrary scale A. 

An important non-trivial test of our framework is to 
calculate the /3-function of the WZ model which is well 
known to four loop order. We start with equation 0. 
Because go cannot depend upon our arbitrary scale A we 
are led to 



35 dZ^ _ 3g ■ 



(54) 



where the dot stands for the logarithmic derivative, i.e. 
A = X{d/dX)A. Now we revoke (|10|l to write a coupling 
constant expansion for the /3 function, viz. 



'Mir 



(Z3 + 5Z1Z2 + 2Z1Z2+4Z: 

+ • • • . (55) 

Using our results displayed in equations H2U|) . (|30|) and 
(ISSll and that 



r(3) 



-2& , 

2Iiog{>^'^) - 26 and 

(2), 



log ^ 



(56) 



we get 




(57) 



Some comments are in order. Whilst the one and 
two-loop coefficients of the f3 function are scheme- 
independent as opposed to the three loop coefficient, the 
coefficient of ^(3) should be universal. That is because 
the diagram from which such coefficient stems is prim- 
itive. Our results for the /3-function exactly agree with 
the calculation performed within differential renormal- 
ization 29]. In contrast, the coefficient we obtain for 
C(3) is different from the one obtained in dimensional 
regularization in p5| . Later, an extension of the work of 
ref. |25j to four loop order in dimensional regularization 
1^ calculated the same coefficient of ^(3) as the one ob- 
tained by us. They have also verified that their result was 
in agreement with a consistency condition which relates 
the coeficients of Z4 and Z3. Finally the results obtained 
by us here and in differential renormalization ^29j for the 
three loop (scheme-dependent) coefficient corresponds to 
a momentum space subtraction scheme (MOM) in which 
a subtraction is perfomed at = A^ |2^] . 



VIII. CONCLUSIONS AND PERSPECTIVES 

If on one hand at the near future experiments at the 
LHC or at a linear e+e~ collider will test decisively su- 
persymmetric extensions of the standard model, on the 
other hand the theoretical machinery which exploits such 
predictions must have a thorough control upon the regu- 
larization and renormalization of supersymmetric Yang- 
Mills theories. 

Since complete regularization framework which is both 
gauge and supersymmetric invariant has yet not been 
constructed, it is reasonable to exploit a framework which 
is gauge invariant and works in the physical dimension 
of the theory. This work is a fundamental step in this 
direction which should help to address important issues 
such as; 

1. gauge field theories with soft supersymmetry- 
breaking terms: for instance QCD with soft break- 
ing which has a particular phenomenological inter- 
est. In jsj models with soft-breaking terms have 
been studied in the Wess-Zumino gauge. However 
there appears new parameters which have no clear 
interpretation as being either susy or soft breaking 
terms (see however U|] ) . IR should be an useful tool 
to renormalize softly broken susy gauge theories 
avoiding some complications used in dimensional 
methods such as the introduction of Slavnov- Taylor 
identities as constraint equations in order to ensure 
both gauge invariance and supersymmetry [3^ . 

2. the anomaly puzzle: the axial and energy momen- 
tum trace anomalies seem to lead to the conclusion 
that the P function of supersymmetric gauge field 
theories should be exhausted by the first loop cor- 
rection Whilst this appear to be the case in 
models with N = 2 supersymmetry, the case A'^ = 1 
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is less clear 34] particularly when one uses dimen- 
sional regularization methods. In D£R which 
like IR does not resort to analytic continuation on 
the space time dimension, the infrared origin of the 
two loop coefficient of the /^-function of = 1 SYM 
has been discussed (see also [2^). IR, which oper- 
ates in momentum space and on the physical di- 
mension of the theory, clearly isolates infrared and 
ultraviolet divergences by means of distinct scales. 
Thus we expect that IR can shed some light on 
the infrared effects and scheme dependence of the 
higher order corrections (if any) to the /3-function 
of TV = 1 supersymmetric gauge theories as well as 
establish a correspondence with other frameworks 
such as the canonical Wilsonian coupling flow . 



discussed earlier it should be infrared safe. To see that 
explicitly let us split the integral from to oo in three 
pieces: from to 7^ plus from to p^ (infrared 
cutoff) plus a remaining contribution from /i^ to 00. The 
latter can be easily shown to be related to H28I) . 



fi '■>(—) 

J^2 X \ X J 



1 



(64) 



whereas the two other pieces added together yield 
l/2{ln2[pV(A2e2)] - / [X^e')]} . Note that ^such 

contribution diverges as /i — > 0. However this behaviour 
is tamed if we correctly define a genuine ultraviolet di- 
vergent object by exchanging /i against A 7^ with the 
help of equation 129|) . Putting all the results together we 
arrive at 



Appendix A 

We can easily obtain the finite part of G{p^) in a simple 
form using Rosner's technique of expanding propagators 
in Chebyshev's polynomials [? ]. 

In Euclidean space H25|) reads: 



-(4^) 



1 , /A2e2 
In 



1E 



(58) 



Because Ge{p^) is a function of p^ only we may average 
over the directions of p in four dimensions. Let C,„ be 
the Chebyshev's polynomial of order m. It satisfies a four 
dimensional version of its orthogonality relation, viz. 



dVL 
2^2 



^ V BO / \ VQ / 



pq / \ pq 
with Co — I- Following [? ] we write 



1 



{p + qf 



pq ^-^ \ pq J 



where 



p/q iip<q 



q/p ilp> q 
Using ^ and ^ in ^ yields 

1 



Ge{p') = -{AiiY — (<zb)lnl 2 

qE I P ^ T 



A2e2 



(59) 



(60) 



(61) 



(62) 



/^^ = 1/(271)^ J dQq J q'^ dq. After integrating over the 
angles we get 



GB(p2) = -ln 



AV 



- 1- 



-In 



dx. (63) 



Take the last term in the equation above. That integral 
contains the ultraviolet divergent piece from which we 
define the renormalization constant Z2. Moreover as we 



G.(p^) = -ln(^)-ihr 



p2 J b 



, (65) 



from which we may finally define the finite part of the 
two loop amplitude which in turn satisfies a Callan- 
Symanzik renormalization group equation with renormal- 
ization scale A. In Minkowski space it reads: 



Gfi„(p2) = -ln 



\ 2 2 

-A e 



1, 2/-AV 
- m 



(66) 



Appendix B 

The graph displayed in figure El contains no subdia- 
grams. The irreducible loop integral that represents the 
ultraviolet divergent content of this graph can be dis- 
played by studying in equation . It will be 
convenient to express the corresponding counterterm as 
one of the basic divergent integrals which have appeared 

so far, namely Iiog[\^), Iioli^'^), etc. Because 

this diagram is primitive and ultraviolet logarithmically 
divergent, is expected to be proportional do Iiog{X^)- 
We use the same technique of appendix A in order to 
extract such ultraviolet behaviour. Moreover, it will be- 
come clear why the coefficient of ^(3) in (|52|) . which will 
appear in the three loop contribution to the /3 function, 
is universal (save a coupling constant redefinition involv- 
ing C(3)) despite the /^-function being scheme dependent 
at the three loop level. 

We temporarily work in Euclidean space where 
reads 



H, 



2 A,,fc9W(/-g)2(/-fc)2(g_A:)2 
Equations H7U|) and (|60|l enable us to write 



(67) 



(68) 
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where 



y ^ i?h [ E {i\qr+'{i\kr+'{Q\kK' 



27r2 J 27r2 g fc 

m,n,p 

X Cnil,q)Cmil,k)Cp{q,k) 



(69) 



Thus the identity 

appUed to yields, after straightforward algebra, 

n=0 ' 

Therefore, in Minkowki space, after introducing an in- 
frared cutoff /i and making use of equation (|18(l we obtain 



= -3z(4^)^C(3)/io3(A^) + 3C(3) In - , (72) 



where lim^^o[3C(3) ln(AV/i2) + ^6(p', m')] = Geip^X') 
is well-defined. 



Acknowledgements 

The authors wish to thank Prof. Victor Rivelles for en- 
lightening discussions and CNPq-Brazil for the financial 
support. M.C.N, thanks S. Gobira for useful discussions 
during early stages of this work. 



[1] N. N. Bogolyubov, D. V. Shirkov, Introduction to the 

theory of the quantised fields, Wiley, 1980; K. Hepp, [16 
Comm. Math. Phys.2 (1966) 301; W. Zimmermann, Lec- 
tures on Elementary Particles and Quantum Field The- [17 
ory, Brandeis (1970), edited by S. Deser, M. Grisaru and 
M. Pendleton; J. H. Lowenstein, BPHZ Renormalisa- [18 
tion, Proceedings of NATO Advanced Study Institute, 
edited by G. Velo and A. S. Wightman, Erice, 1975. [19 
J. H. Lowenstein and B. Schroer, Phys. Rev. D6 (1972) 
1553; O. Piguet and K. Sibold, Nucl. Phys. B248 (1984) [20 
301. 

W. Siegel, Phys. Lett. B84 (1979) 193; idem, Phys. Lett. [21 
B94 (1980) 37. 

W. Hollik and D. Stockinger, Eur. Phys. J. C20 (2001) 
105. [22 
See for instance: I. Jack and D. T. R. Jones in "Perspec- 
tives in Supersymmetry" , Editor G. Kane, World Scien- [23 
tific, and references therein. 

D. Z. Freedman, K. Johnson and J. I. Latorre, Nucl. [24 
Phys. B371 (1992) 353. 

P. E. Haagensen and J. I. Latorre, Ann. Phys. 221 (1993) 
77. [25 
M. Chaichian and W. F. Chen, Phys. Lett. B409 (1997) 
325. [26 
W. F. Chen, H. C. Lee and Z. Y. Zhu, Phys. Rev. D55 
(1997) 3664. [27 
D. Z. Freedman, G. Grignani, K. Johnson and N. Rius, 
Ann. Phys. 218 (1992) 75. 

F. del Aguila, A. Culatti, R. Munoz Tapia and M. Perez- [28 
Victoria, Nucl. Phys. B537 (1999) 561; Nucl. Phys. [29 
B504 (1997) 532 and hep-th/9711474 in Barcelona 1997, [30 
Quantum effects in the minimal supersymmetric stan- 
dard model, 361-373. [31 
M. Perez- Victoria, Phys. Lett. B442 (1998) 315. 
F. del Aguila and M. Perez- Victoria, Acta Phys. Polon. [32 
B29 (1998)2857. 

P. E. Haagensen and J. L Latorre, Phys. Lett. B283 [33 
(1992) 293; idem Phys. Lett. B419 (1998) 263. 
F. del Aguila and M. Perez- Victoria, |hep-ph/990'T29Tl [^^ 



and Acta Phys. Polon. B28 (1997) 2279. 

T. Hahn and M. Perez- Victoria, Comput. Phys. Com- 

mun. 118 (1999) 153. 

O. A. Battistel, PhD thesis. Federal University of Minas 
Gerais, Brazil (2000). 

O. A. Battistel, A. L. Mota and M. C. Nemes, Mod. Phys. 
Lett. A13 (1998) 1597. 

A. P. Baeta Scarpelli, M. Sampaio and M. C. Nemes, 
Phys. Rev. D63 (2001) 046004. 

A. P. Baeta Scarpelli, M. Sampaio and M. C. Nemes, B. 

Hiller Phys. Rev. D64 (2001) 046013. 

M. Sampaio, A. P. Baeta Scarpelli, B. Hiller, A. Brizola, 

M.C. Nemes, S. Gobira, Phys. Rev. DD65 (2002) 

125023. 

A. Brizola, O. A. Battistel, M. Sampaio and M. C. 
Nemes, Mod. Phys. Lett. A14 (1999) 1509. 

0. A. Battistel and M. C. Nemes, Phys. Rev. D59 (1999) 
055010. 

S. R. Gobira and M. C. Nemes, Perturbative n-loop renor- 
malization by an implicit regularization technique, to ap- 
pear in Int. J. of Mod. Phys. (2003). 
L. F. Abbott and M. T. Grisaru, Nucl. Phys. B169 
(1980) 415. 

M. T. Grisaru, M. Rocek and W. Siegel, Nucl. Phys. 
B159 (1979) 429. 

Prem P. Srivastava, Supersymmetry, Superfields and Su- 
pergravity, an introduction. Graduate Studies Series in 
Physics (1986). 

G. Dunne, N. Rius, Phys. Lett. B293 (1992) 367. 
P. E. Haagensen, Mod. Phys. Lett. A7 10 (1992) 893. 
A. Sen and M. K. Sundaresan, Phys. Lett. BlOl (1981) 
61. 

N. Maggiore, O. Piguet and S. Wolf, Nucl Phys. B476 
(1996) 476. 

M. Sampaio, M. C. Nemes, A. P. Baeta Scarpelli , Work 
in Progress. 

V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. 

1. Zakharov, Phys. Lett. B157 (1985) 169. 

L Jack, D. Jones and C. North, Nucl. Phys. B473 (1996) 
308; O. V. Tarasov, V. A. Vladimirov, Phys. Lett. B96 



10 



(1980) 94; M. T. Grisaru, M. Rocek and W. Siegel, Phys. 
Rev. Lett. 45 (1980) 1063; W. Caswell and D. Zanon, 
Phys. Lett. BlOO (1980) 152. 
[35] J. Mas, M. Perez- Victoria and C. Seijas, JHEP 0203 
(2002)049 



[36] A. A. Soloshenko and K. V. Stepanyantz, Three loop P 
function for N = 1 supersymmetric electrodynamics reg- 
ularized by higher derivatives, liep-th/0304083 

[37] N. Arkani-Hamed and H. Murayama, JHEP 0006 (2000) 
030. 



